Mutual adaptation of a Faraday instability pattern with its flexible boundaries in 

floating fluid drops* 
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Hydrodynamic instabilities are usually investigated in confined geometries where the resulting 
spatio-temporal pattern is constrained by the boundary conditions. Here we study the Faraday 
instability in domains with flexible boundaries. This is implemented by triggering this instability in 
a floating fluid drop. An interaction of Faraday waves with the shape of the drop is observed, the 
radiation stress of the waves exerting a force on the surface tension held boundaries. Two regimes 
are observed. In the first one there is a co-adaptation of the wave pattern with the shape of the 
domain so that a steady configuration is reached. In the second one the radiation stress dominates 
and no steady regime is reached. The drop stretches and ultimately breaks into smaller domains 
that have a complex dynamics including spontaneous propagation. 

PACS numbers: 05 65. +b, Self-organized systems, 05 45. -a, Non linear dynamics and chaos, 47. 20. Ma 
Interfacial instabilities , 47.35.Pq, Capillary waves 



Fluid dynamics instabilities usually appear in two 
types of situations corresponding to confined or opened 
geometries respectively. For instance thermal buoyancy, 
when confined in a box, gives rise to patterns of Rayleigh- 
Benard rolls adapted to the boundaries pQ . In contrast in 
an open medium, an isolated source of heat generates a 
thermal plume in which the growing turbulent structures 
define the envelop of the unstable region [2 , 3 . We ad- 
dress an intermediate situation in which an instability de- 
velops inside a finite domain with flexible boundaries and 
we study the interplay between the pattern and its bor- 
ders. We use the Faraday instability in which waves form 
on the surface of a fluid submitted to vertical oscillations 
[4] and whose resulting patterns have been widely studied 
in confined geometries [SJ |6] . Here we choose to confine 
the Faraday instability in a drop of low viscosity float- 
ing on a very viscous, stable, immiscible fluid [Figjlja)]. 
We use a classical Faraday experiment set-up. A circular 
cell (radius 10cm, depth 0.8cm) is placed on a vibration 
exciter generating a vertical oscillation of acceleration 
j(t) = 7 m cos(2tt fat). The investigated frequency and 
amplitude ranges are respectively 50Hz < fo < 250Hz 
and < jm/g < 10, where g is the acceleration of grav- 
ity. The motion can be observed by a stroboscopic video 
camera or a high speed video camera. The heavier and 
more viscous fluid (fluid 1) fills the cell and forms a bath 
of thickness 5mm. A controlled quantity V2 of the less 
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FIG. 1: (a) Vertical section of an isopropanol drop floating 
on perfluorated oil in absence of oscillations, (b-d) A drop of 
volume V 2 = 1.00 ± 0.02ml forced at f = 130Hz for three 
amplitudes of forcing, (b) Faraday waves of weak amplitude, 
Jm/g = 3.28, (c) steady elongated state, jm/g = 3.93 and (d) 
highly perturbed state jm/g = 6.48. (e) A drop of volume 
V2 = 10.0 =b 0.1ml square when forced at fo = 160Hz with 
7m/ ' g — 4.90. The bars are 1cm long. 
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viscous fluid (fluid 2) is then deposited becoming a float- 
ing drop that forms a single circular pancake far from the 
boundaries. Its radius ro and thickness Y12 at rest result 
from the equilibrium between buoyancy and capillarity 
P3 [8] . The size of the bordering meniscus fixed by the 
capillary length is small compared to the horizontal size. 
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In all our experiments we find ft 2 — 2mm and Tq ~ 1.3cm. 

There is a range of forcing amplitude for which the 
Faraday instability forms in the drop only. The resulting 
waves are initially disordered and generate fluctuations of 
the drop boundary [Figfljb)]. In return the fluctuations 
of the boundary generate an unsteadiness of the wave 
pattern. Two possible archetypes of evolution can then 
be observed [Fig{l][2] that depend on the fluids. We in- 
vestigated many pairs of immiscible fluids, having differ- 
ent viscosity contrasts and wetting properties, and found 
only these two possible behaviours. 

In the first archetype, the system self-organizes by a 
co-evolution of the wave field and the boundaries so that 
an equilbrium is reached [9] . The example shown in figjl] 
was obtained with isopropanol (density p2 = 785kg/m , 
viscosity P2 = 1.8 x 10 _3 Pa • s and surface tension against 
vapor (J2 = 24mN/m) floating on perfTuorated oil (pi = 
1850kg/m 3 , pi = 26 x 10" 3 Pa • s and a x = 15mN/m). 
The measured interfacial tension is = 6.3mN/m. 
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FIG. 2: The case of an ethanol drop of volume V{ — 1.00 =b 
0.02ml deposited on silicon oil. (a) Sketch of the vertical 
section of the floating drop in the absence of oscillations, (b) 
It is circular at rest with area and thickness A2 = 459 =b 5mm 2 
and ft 2 = 2.2mm respectively, (c-e) Three successive images 
showing the temporal evolution of the drop when forced at 
fo — 130Hz with jm/g = 7.00. The bar is 1cm long. 

In the second type of evolution the displacement of the 
boundary due to the waves does not lead to equilibrium 
[TO] . The formation of parallel standing waves results 
into a constantly increasing elongation of the drop into a 
snake-like structure which breaks into fragments having a 
large variety of dynamical behaviors. The example shown 
in fig (2] was obtained with ethanol (p' 2 = 789kg/m 3 , 
p' 2 = 0.9 x 10" 3 Pa-s, cr' 2 = 23mN/m) floating on sil- 
icon oil (p[ = 965kg/m 3 , p[ = 100 x 10 _3 Pa • s and 
a[ = 20mN/m). The measured interfacial tension is 



cr' 12 0.7mN/m. 

The difference between the two regimes can be under- 
stood by a dimensional analysis evaluating the ratio of 
the destabilizing factor (the radiation pressure) to the 
restoring one (capillary pressure) as will be discussed be- 
low. 

We first characterize experimentally the former 
archetype as a function of the frequency and the forc- 
ing amplitude. A phase diagram depicting the different 
states of the system is shown in figj3] When the forc- 
ing amplitude is increased, Faraday waves appear in the 
drop above a first onset (j m > 7^). They form complex 
unsteady patterns which deform the drop [FigQb)], its 
average shape remaining circular. Above a second thresh- 
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FIG. 3: Phase diagram of the different regimes observed in 
the case of an isopropanol drop of volume V2 = 1.00 ±0.02 ml 
deposited on perfluorated oil. The drop area at rest is A 2 = 
531 =b 11mm 2 , its average thickness ft 2 — 1.9mm. C = circular, 
D = deformed, E = elongated, HP = highly perturbed, F = 
Faraday instability in oil. 

old 7^ the waves strengthen and start organizing, result- 
ing into an elongation of the drop in the perpendicular 
direction. Correlatively this elongation favors a further 
organization of the waves. By a slow evolution the drop 
reaches a stable elongated shape and the waves become 
steady too [Figjljc)]. They are perpendicular to the long 
sides of the drop (as usual for Faraday waves in elongated 
cells). By increasing the forcing the drop undergoes a fur- 
ther elongation to a new steady state. The shape can be 
defined by the ratio R of the length of the minor over 
the major axis of the drop. FigQa) shows the evolution 
of R as a function of j m for various frequencies. Above 
a third threshold 7^ p , the waves become chaotic and 
the envelope of the domain undergoes very large fluctua- 
tions [Figjljd)]. Finally at a fourth threshold 7^ Faraday 
waves form on the viscous substrate. We cannot obtain 
the second archetype by increasing the forcing amplitude. 
We measured the area A2 covered by the drop during this 
process and found that it does not change. The observed 
wavelength is given by the dispersion relation of gravity- 
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capillary waves for the Faraday frequency fo/2. 

The first type of behavior can be understood as result- 
ing from the effect of the radiation pressure P r [TT] of the 
surface waves distorting the boundaries of the pancake. 
We thus seek stable solutions in which the radiation and 
hydrostatic pressures are balanced by capillarity. We as- 
sume an unidirectional standing wave along the x-axis, 
with small thickness for the drop and small wavelength 
for the wave, the initial radius of the drop being the 
length unit. The small thickness assumption allows a lu- 
brication approximation transforming the 3D boundary 
value problem in a 2D one while the small wavelength 
approximation allows an average on larger length-scales. 
An asymptotic analysis adapted from [12] shows that the 
function y(x) describing the drop shape satisfies a two- 
dimensional Laplace law modified by the radiation pres- 
sure along the normal. It is a Riccati equation: 
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where Ph is a term of hydrostatic pressure (assumed con- 
stant for quasi-steady boundary), f{0) depends on the 
wetting angle and P r is the radiation pressure 
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with A the amplitude of the Faraday waves, uo — tt/q 
the angular frequency and F = e hl ~ h2 (1 + B cosh — 
Bsmhkhi, with B = —k/uu 2 [(p2 — pi)/p2g — &i,2k 2 /pi) 1 a 
factor that takes into account thickness effects and waves 
transmission for a baroclinic mode in fluid 1, k being the 
wave vector modulus. Using the volume conservation, 
the parameters of eq{l] reduce only to one free parameter 
responsible for the final equilibrium shape, that is the 
ratio a = P r /Ph. In this case we have found explicit 
analytical solutions of eq(T] The drop shape is given by 
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sin 2 $0 — sin 2 P^y/ax 



with sin 2 $ = a/(l + a) and P' h = P h /[a 2 f(0)}. The con- 
servation of the area gives P' h = [log(l + a) /ay/1 + «] 1//2 , 
being P^ = 1 at rest (P r = 0). When R > 0.25 good 
fits of the experimental shapes are obtained [FigQb, c)]. 
The evolution of the aspect ratio R with forcing can also 
be obtained analytically 
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We measured the amplitude A of the waves in the drops 
by lateral zoomed-in films and checked that it varies as 
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m) 1 / 2 - Using it to evaluate P r and thus a, we 



FIG. 4: (a) Evolution of the aspect ratio of an elongated 
isopropanol drop on perfluorated oil as a function of the forc- 
ing amplitude for different frequencies. The lines are the fits 
by eq|4] (b, c) Photographs of two drops in an elongated 
state at f = 130Hz for j m /g = 3.97 and R = 0.42 (b) and 
lm/g = 4.62 and R — 0.28 (c) and the predicted shapes by 
the corresponding solutions [Eq|3]. The bar is 1cm long. 



obtain fits of i?(7 m ) for various frequencies [FigQa)]. 
The computed shapes become inconsistent with the ini- 
tial small wavelength and small thickness hypothesis of 
the model whenever the radius of curvature at the tip 
of the computed shape becomes small. For this reason 
the very elongated shapes (R < 0.25) obtained with this 
model cannot account for reality. Experimentally we ob- 
serve that the shape of the tip becomes fixed while the 
drop keeps elongating, the curvature radius of the tip 
being of the order of the wavelength. 

We can now turn to the second archetype. The main 
characteristic is that ethanol is wetted by silicon oil. At 
rest an oil film is observed to cover the upper surface of 
the drop [Figj2ja)]. The instability in ethanol appears 
in a subcritical way i.e. large amplitude waves form 
at threshold. These waves stretch the drop [Figj2^c)] 
but there is no convergence towards a final stable shape. 
There is no theory for this dynamical regime yet. 

However its existence can be understood by dimen- 
sional analysis. We use ao = P r /P®, the ratio of the 
estimated wave radiation pressure [Eq{2] to the 2D pres- 
sure due to capillary effects in the circular drop at rest. 
For drops of centimetric size and Faraday waves of the 
usually observed amplitude this analysis applied to the 
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first archetype gives ao — 0.1. In the second archetype, 
because of the wetting, the capillary tension is the inter- 
facial one <j[ 2 = 0.7mN/m. As a result we find a' ~ 2. 
While in the first case the pressure exerted by the waves 
is only a perturbation of the capillary equilibrium of the 
drop, in the second this radiation pressure exceeds the 
possible response of capillary forces so that no steady 
solution can be reached. 

The elongation is followed by buckling [FigMd)] then 
breaking into several fragments [Fig{2|e), Figj5ja)]. We 
measured the total area A' 2 covered by the pattern and 
found that it increases with forcing, proving that the 
drops are stretched by the waves (their thickness de- 
creases). By buckling large fragments can take croissant 
or a horseshoe shapes [FigjEJd)]. These shapes are sta- 
tionary but propagate in the direction of the curvature. 
This is one more example of self-propagation by a spon- 
taneous symmetry breaking [13] . The velocity is constant 
and a croissant (or a horseshoe) could move indefinitely 
in an infinite bath. In practice a bath of finite size is 
covered with moving and motionless fragments [FigjEJa)] 
that keep colliding, merging and splitting. The global 
aspect is reminiscent of the interplay of structures ob- 
tained in cellular automata such as Conway's game of 
life [14]. Very small fragments are observed to remain 
steady with a stationary shape that recalls the elongated 
one of the first regime [FigJEJb)]. This possibility of a 
steady regime for small drops can be understood by the 
dimensional analysis : ao becomes smaller when the ra- 
dius of the drop at rest is very small. 
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FIG. 5: (a) The surface of the experimental cell after the frag- 
mentation of a drop of ethanol on silicon oil for fo — 80Hz and 
Jm/g = 5.06. (b-d) The shape of some of the resulting steady 
fragments. Each one corresponds to a different fragment size. 
The bars are 1cm long. 

A remarkable feature of the snake-like structure is that 
they have constant transverse widths and tip radii, both 
being of the order of the Faraday wavelength. Their 
stretching and instabilities affect mostly the middle of 
their length. The buckling could be due to the generation 



of a streaming flow [15] that has been already observed in 
Faraday instability [16] . This streaming effect could also 
be responsible for their motion on the substrate. The in- 
vestigation of the resulting dynamical regimes is beyond 
the scope of this letter. 

We have found that, by slow dynamics, a mutual adap- 
tation is possible between an instability and its bound- 
aries. This phenomenon is related to the self-tuning of 
oscillators. In optical cavities the radiation pressure was 
shown [17] to create a coupling between the mirrors de- 
grees of freedom and the optical field. Similarly, a forced 
mechanical oscillator with an additional degree of free- 
dom [18 exhibits a slow drift by which it can self-tune. 
The present phenomenon is more general and should 
show up in other type of instabilities confined in adapt- 
able boundaries. 
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